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Critical finite size scaling functions for the order parameter distribution of 
the two and three dimensional Ising model are investigated. Within a recently 
introduced classification theory of phase transitions the universal part of crit- 
ical finite size scaling functions has been derived by employing a scaling limit 
which differs from the traditional finite size scaling limit. In this paper the 
analytical predictions are compared with Monte Carlo simulation results. We 
find good agreement between the analytical expression and the simulation res- 
ults. The agreement is consistent with the possibility that the functional form 
of the critical finite size scaling function for the order parameter distribution 
is determined uniquely by only a few universal parameters, most notably the 
equation of state exponent. 
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A universality class in the theory of critical phenomena is identified by a set of critical 
exponents and a set of universal scaling functions [1-3]. In practice finite size scaling functions 
and their associated universal amplitudes or amplitude ratios have become an indispensable 
tool for the extraction of universal behaviour from numerical simulations of finite systems 
throughout many fields of physics [4-8], and the subject remains of vigorous research interest 
[9-12]. It is therefore of broad interest to obtain exact information on finite size scaling 
functions. Such knowledge would also be particularly important for discerning corrections 
to scaling behaviour. 

Most determinations of critical finite size scaling functions or amplitude ratios have been 
carried out by numerical simulation of critical systems [3,13-21]. Analytical calculations 
[22-25] are complicated because at criticality the influence of boundary conditions cannot be 
neglected. This renders renormalization group methods difficult to implement. The problem 
of hyperscaling violations [26,27] has further obscured the basic question as to what extent 
finite size scaling functions are universal or not. Recently [28], however, finite size scaling 
theory has been reanalysed from the perspective of a general classification theory for phase 
transitions [29-33]. In these papers it was shown that phase transitions in statistical mech- 
anics may be classified according to a generalized classification theory. This theory classifies 
each transition according to its generalized order in analogy with Ehrenfest's classification 
scheme. The objective of the present letter is to compare analytical predictions of the gener- 
alized classification scheme with computer simulation results. The comparison is carried out 
for the critical finite size scaling function of the order parameter distribution in the two- and 
three dimensional Ising models with periodic boundary conditions. Good agreement between 
simulations and the theoretical predictions is found over the entire range for which numerical 
data are available. 

Define p(ip, L } £) to be the probability density function for the fluctuating order parameter 
^ in a finite system of size L and order parameter correlation length £. Then the scaling 
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function of interest p(x, y) is defined by 



p(ll>,L,{) = l d ( d t-d")/(d-d t ) p^L^-nKd-n^l^) 



(1) 



where dq, is the anomalous or scaling dimension of the order parameter, d* is Fishers an- 
omalous dimension of the vacuum [34], and is Binders thermodynamic length [35]. If 
hyperscaling holds then d* = 0, the thermodynamic length becomes the correlation length, 
£o = £ 5 an d the exponent in (1) reduces to the familiar form dq, = (3/u where (3 is the order 
parameter exponent and v the correlation length exponent. The scaling function p(x } y) is 
expected to be universal up to the choice of boundary conditions [13]. 

Given the scaling Ansatz (1) the traditional scaling analysis [13,25] of the critical scaling 
function p(x,0) distinguishes two cases, (i) For i < 1 the scaling function is expected to 
have the universal Landau-Ginzburg form [13] 



(ii) For x ^> 1 and groundstate boundary conditions (e.g. all spins positive for the Ising 
model) the scaling function is expected to have the squeezed exponential form [36] 



where 8 denotes the equation of state exponent. Based on scaling arguments the same form 
is expected to apply for periodic boundary conditions [25]. 

Little exact information is available for p(x, 0). To the best of our knowledge only the cu- 
mulant ratio g(0) = (/ \x\ 4 p(x } 0) dx) / (f \x\ 2 p(x } 0) dx) 2 for the two dimensional Ising model 
with singular periodic boundary conditions has been calculated exactly [24]. For noncritical 
systems on the other hand the analogous noncritical scaling function is Gaussian by virtue of 
the central limit theorem [13]. The absence of exact information about the critical function 






(3) 
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p(x,0) even for the otherwise exactly solvable two dimensional Ising model is related to the 
absence of its solution in nonzero magnetic field [36]. 

Recently the universal part of critical finite size scaling functions has been related to finite 
ensemble scaling functions [33,28]. Finite ensemble scaling functions arise in the ensemble 
limit while finite size scaling functions arise in the finite size scaling limit. For a (/-dimensional 
discretized lattice system in the fully finite hypercubic geometry the finite size scaling limit is 
defined as the limit z,, 1 ™^ in which the box dimension L and the correlation length £ increase 

L/( = c 

to infinity in such a way that their ratio remains constant. In the ensemble limit on the other 
hand the lattice constant a approaches simultaneously. More precisely, the ensemble limit 
is defined as the limit m,n™oo in which M = (£/ ' a) d and N = (L/£) d approach infinity such 

N/M=c 

that their ratio remains constant. 

In the finite size scaling limit the critical finite size scaling functions are found to contain 
a universal as well as a nonuniversal part. The universal part is given by the finite ensemble 
scaling functions which arise in the finite ensemble limit and for which analytical expressions 
can be derived if the scaling dimension of the critical operator in question is known. For a 
critical operator X at the critical point of a (/-dimensional system with Ising symmetry and 
periodic boundary conditions the universal part of the critical finite size scaling function is 
written as [28] 

p(x,0) = ^h + (x;w x ) + ^h~(x;w x ) (4) 

with 

wx = min(2, 2 — ax) (5) 

where ax is the thermodynamic fluctuation exponent [34] of the observable X. If X rep- 
resents the energy density then as = a, the specific heat exponent, while for the order 
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parameter density ctq, = 1 — (l/S) where 8 is the equation of state exponent. The scaling 
functions h ± (x; wx) obey h + (x) = h~(—x) and can be written in terms of the ^-function 
representation of stable probability densities [37-39] as 



The general class of _ff-functions is usually defined in terms of Mellin-Barnes contour integrals 
and contains Meijer's G-function as well as many other generalized hypergeometric functions 
as special cases. For a precise definition we refer the reader to standard tables [38]. Note that 
equation (6) depends only on wx which is completely determined by the scaling dimension 
of X. Note also that equations (4) and (6) apply for periodic boundary conditions. The 
general theory identifies a universal shape and symmetry parameter which is related to 
different choices of boundary conditions [28]. It was identified in [28] to be unity for periodic 
boundary conditions and the same choice has been applied here. In the following we focus on 
order parameter fluctuations, i.e. X = ^ . In that case the index wx becomes wq, = 1 + (1/8) 
where 8 is the equation of state exponent. 

The scaling function given by (4) and (6) are consistent with the scaling results (2) and 
(3). The functions h ± (x;wx) are entire functions of x, and thus p(x,0) may be expanded 
around x = as assumed in (2). Secondly the asymptotic expansion of the _ff-functions [40] 



which is consistent with the scaling result (3) even though it is not derived in the same 
scaling limit [41]. 

In Figure la we compare the analytical result to simulation data for the two dimensional 
Ising model with periodic boundary conditions. In this case 8 = 15 and the critical tern- 
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perature is known exactly. The simulation results are represented as crosses, the analytical 
result as a solid line. The crosses give a smoothed representation of scaled simulation data 
[19] for system size L = 64. All distributions are scaled to unit norm and variance. For 
the analytical curve this requires a cutoff which was chosen at a value close to the largest 
simulation data point. 

In Figure lb the same comparison is shown for the case of the d = 3 Ising model. In 
this case neither the critical temperature nor the equation of state exponent 8 are known 
exactly. The data points represent original high precision Monte Carlo simulations in which 
systems of size L = 20 and L = 32 were studied for 50 X 10 6 Monte Carlo sweeps each using 
a vectorized code on a Cray YMP. The data for the system size L = 20 is represented by 
crosses, that for L = 32 by circles. We used the estimates [20] J/(kBT c ) = 0.2216595 for the 
critical temperature and 8 = 4.8 for the equation of state exponent. Here J denotes the Ising 
exchange coupling and ks is Boltzmann's constant. 

The agreement between the analytical prediction and the simulation results in both two 
and three dimensions is gratifying. We attribute the small discrepancies, in part at least, to 
the low statistics in the tails of the scaling function. This view is supported by comparing the 
scaling function obtained in a previous small-scale study of the 3D Ising model [13], with that 
of figure 2. While both the scaling functions of [13] and those reported here exhibit excellent 
data collapse, the scaling functions in both cases are markedly different, demonstrating the 
presence of the nonuniversal part [28] . The good agreement between theory and simulations 
lends substantial support to the theoretical ideas from which the scaling functions derive. In 
particular it should be emphasized that the only difference between the solid lines in Figure 
la and lb is the value of 8. This suggests that the full functional form of the universal 
part of the critical finite size scaling function of the order parameter distribution could in 
simple cases be determined by a few universal parameters, most notably the equation of state 
exponent. Further numerical and analytical studies are, however, required to conclusively 
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establish whether the proposed universal scaling functions agree fortuitously with the Monte 
Carlo data or are more generally correct. 
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FIGURE CAPTIONS 



FIG 1 : Comparison of the analytical predictions of equations (4)-(6) with Monte-Carlo 
simulation data for the magnetisation distribution of the two and three dimensional critical 
Ising model with periodic boundary conditions. In (a) the 2D data of [19] (crosses) is 
compared with the analytical prediction for 8 = 15 (full curve). Part (b) shows original 
3D simulation data for L = 20 (crosses) and L = 30 (circles), collected at temperature 
J I '(ksTc) = 0.2216595. The corresponding analytical curve corresponds to the estimate 
8 = 4.8 [20]. All data have been scaled to unit norm and variance and statistical errors do 
not exceed the symbol sizes. 
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